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We use recent observations of the Hi-mass function to constrain galaxy formation. 
The data conflicts with the standard model where most of the gas in a low-mass dark 
matter halo is assumed to settle into a disk of cold gas that is depleted by star for- 
mation and supernova-driven outflows until the disk becomes gravitationally stable. 
. Assuming a star formation threshold density supported by both theory and observa- 

' tions this model predicts HI masses that are much too large. The reason is simple: 

supernova feedback requires star formation, which in turn requires a high surface den- 
1/"") 1 sity for the gas. Heating by the UV background can reduce the amount of cold gas in 

^sO ' haloes with masses < 1O 9 5 /i _1 M0, but is insufficient to explain the observed HI mass 

function. A consistent model can be found if low-mass haloes are embedded in a pre- 

■ heated medium, with a specific gas entropy ~ lOkeVcm 2 . In addition, such a model 

simultaneously matches the faint-end slope of the galaxy luminosity function without 

. the need for any supernova driven outflows. We propose a preheating model where 

(*S the medium around low-mass haloes is preheated by gravitational pancaking. Since 

gravitational tidal fields suppress the formation of low-mass haloes while promoting 
that of pancakes, the formation of massive pancakes precedes that of the low-mass 
haloes within them. We demonstrate that the progenitors of present-day dark matter 
j3 \ haloes with M < lO 12 /^ 1 M Q were embedded in pancakes of masses ~ 5x 10 12 /i _1 M Q 

at z ~ 2. The formation of such pancakes heats the gas to a temperature of 5 x 10 5 K 
and compresses it to an overdensity of ~ 10. Such gas has a cooling time that exceeds 

■ the age of the Universe at z<2, and has a specific entropy of ~ 15keVcm 2 , almost 

exactly the amount required to explain the stellar and HI mass functions. 
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1 INTRODUCTION 

The cold dark matter (CDM) model of structure formation 
has proven a very successful paradigm for understanding 
the large-scale structure of the Universe. However, as far as 
galaxy formation is concerned, a number of important issues 
still remain. A long-standing problem is that CDM models in 
general predict too many low-mass dark matter haloes. The 
mass function of dark matter haloes, n(M), scales with halo 
mass roughly as n(M) oc M~ 2 at the low-mass end. This 
is in strong contrast with the observed luminosity function 
of galaxies, <b(L), which has a rather shallow shape at the 
faint end, with $(!/) oc L" 1 . To reconcile these observations 
with the CDM paradigm, the efficiency of star formation 
must be a strongly nonlinear function of halo mass (e.g. 
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Kauffmann, White & Guiderdoni 1993; Benson et al. 2003; 
Yang et al. 2003; van den Bosch et al. 2003b). One of the 
biggest challenges in galaxy formation is to understand the 
physical origin of this strongly nonlinear relationship. 

Another important but related challenge is to under- 
stand the baryonic mass fraction as a function of halo mass. 
The baryonic mass in dark matter haloes maybe roughly di- 
vided into three components: stars, cold gas, and hot gas. 
In the most naive picture one would expect that each halo 
has a baryonic mass fraction that is close to the univer- 
sal value of ~ 17 percent 1 . In this naive picture one can 
achieve a low ratio of stellar mass to halo mass by keeping 
a relatively large fraction of the baryonic mass hot, either 
by preventing the gas from cooling or by providing a heat- 



1 Corresponding to a ACDM concordance cosmology with 
Q B h 2 = 0.024 and Q m h 2 = 0.14 (Spergel et al. 2003) 
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ing source that turns cold gas into hot gas. Alternatively, 
one may achieve a low star formation efficiency in low-mass 
haloes by making the total baryonic mass fraction lower in 
lower mass haloes, which can be achieved either by blow- 
ing baryons out haloes or by preventing baryons from ever 
becoming bound to haloes in the first place. 

At the present it is unclear which of these scenarios 
dominates. To a large extent this ignorance owes to the 
poor observational constraints on the baryonic inventory as 
function of halo mass. Historically, most observational stud- 
ies of galaxies have focussed on the stellar light. Although 
this has given us a fairly detailed consensus of the relation 
between stellar mass (or light) and halo mass (e.g., Yang 
ct al. 2003, 2005; van den Bosch et al. 2003b; Tinker et 
al. 2004), less information is available regarding the hot and 
cold gas components. Hot, tenuous gas in low mass haloes is 
notoriously difficult to detect, making our knowledge of the 
relation between halo mass and hot gas mass quite limited. 
Based on X-ray observations of a few relatively massive spi- 
ral galaxies, the gas mass in a hot halo component appears 
to be small (e.g. Benson et al. 2000). For cold gas the sit- 
uation has improved substantially in recent years, owing to 
the completion of relatively large, blind 21-cm surveys (e.g. 
Schneider, Spitzak & Rosenberg 1998; Kraan-Korteweg et 
al. 1999; Rosenberg & Schneider 2002; Zwaan et al. 2003, 
2005). Using these surveys, it is now possible to obtain im- 
portant constraints on galaxy formation (see Section 2). 

When modelling galaxy formation, the process most of- 
ten considered to suppress star formation in low mass haloes 
is feedback from supernova explosions. As shown by Dekel 
& Silk (1986) and White & Frenk (1992), the total amount 
of energy released by supernovae can be significantly larger 
than the binding energy of the gas in low mass haloes. There- 
fore, as long as a sufficiently large fraction of this energy can 
be converted into kinetic energy (often termed the 'feed- 
back efficiency'), one can in principle expel large amounts of 
baryonic material from low mass haloes, thus reducing the 
efficiency of star formation. Indeed, semi-analytical models 
for galaxy formation that include a simple model for this 
feedback process are able to reproduce the observed slope 
of the faint-end luminosity function in the standard ACDM 
model, if the feedback efficiency is taken to be sufficiently 
high (e.g., Benson et al. 2003; Kang et al. 2005). 

An important question, however, is whether such high 
efficiencies are realistic. For example, detailed hydrodynami- 
cal simulations by Mac-Low & Ferrara (1999) and Strickland 
& Stevens (2000) show that supernova feedback is far less 
efficient in expelling mass than commonly assumed because 
the onset of Rayleigh- Taylor instabilities severely limits the 
mass loading efficiency of galactic winds. 

This prompted investigations into alternative mecha- 
nisms to lower the star formation efficiency in low mass 
haloes. Another possibility is that photoionization heating 
by the UV background may prevent gas from cooling into 
low-mass haloes (e.g., Efstathiou 1992; Thoul & Weinberg 
1996) by increasing its temperature. Numerical simulations 
have shown that this effect is only efficient in dark matter 
haloes with M < lO lo /i _1 M (e.g. Quinn et al. 1996, Gnedin 
2000; Hoeft et al. 2005), since the gas is only heated to ~ 10 4 
to 10 K. Although this might be sufficient to explain the 
relatively low abundance of satellite galaxies in Milky Way 
sized haloes (Bullock, Kravtsov & Weinberg 2000; Tully et 



al. 2002) , it is insufficient to explain the faint-end slope of the 
galaxy luminosity function. However, if a different mecha- 
nism were to heat the intergalactic medium (hereafter IGM) 
to even higher temperatures (and thus higher entropies), the 
same mechanism could affect more massive haloes as well. 
Such a process is often referred to as 'preheating'. Along 
these lines, Mo & Mao (2002) considered galaxy formation 
in an IGM that was preheated to high entropy by vigor- 
ous energy feedback associated with the formation of stars 
in old ellipticals and bulges and with active galactic nuclei 
(AGN) activity at redshifts of 2 to 3. They showed that such 
a mechanism can produce the entropy excess observed today 
in low-mass clusters of galaxies without destroying the bulk 
of the Ly-a forest. In addition, it would affect the forma- 
tion of galaxies in low mass haloes whose virial temperature 
is lower than that of the preheated IGM. Numerical simu- 
lations show that such preheating may indeed significantly 
lower the gas mass fraction in low mass haloes (e.g., van den 
Bosch, Abel & Hernquist 2003a; Lu et al. in preparation). 

In this paper we investigate an alternative mechanism 
for creating a preheated IGM. Rather than relying on star 
formation or AGN, we consider the possibility that the col- 
lapse of pancakes (also called sheets) and filaments heats the 
gas in these structures and that the low mass haloes within 
them form in a preheated medium. Although the standard 
picture of hierarchical formation is one in which more mas- 
sive structures form later, gravitational tidal fields suppress 
the formation of low-mass haloes, while promoting the for- 
mation of pancakes. Consequently, the formation of massive 
pancakes will precede that of low-mass dark matter haloes, 
which form within them. In this paper we show that the 
shock heating associated with pancake collapse at z < 2 can 
heat the associated gas to sufficiently high entropy that the 
subsequent gas accretion into the low mass haloes that form 
within these pancakes is strongly affected. We demonstrate 
that the impact of this previrialization is strong enough to 
explain both the faint-end of the galaxy luminosity func- 
tion and the low mass end of the galaxy HI mass function, 
without having to rely on unrealistically high efficiencies for 
supernova feedback. 

The outline of the paper is as follows. In § 2, we use 
current observational results of the HI gas mass function to 
constrain star formation and feedback in low-mass haloes. 
We show that these observations are difficult to reconcile 
with the conventional feedback model, but that a consistent 
model can easily be found if low-mass haloes are embedded 
in a preheated medium. In § 3 we describe how shocks asso- 
ciated with the formation of pancakes can preheat the gas 
around low-mass haloes. In § 4 we discuss our results in light 
of existing numerical simulations and discuss the impact of 
our results on the properties of galaxies and the intergalactic 
medium. We summarize our results in § 5. 

2 THE FORMATION OF DISK GALAXIES 
2.1 Observational Constraints 

The models discussed below focus on galaxies that 
form at the centers of relatively low mass haloes with 
M < 10 12 /i _1 M . To constrain these models, we first derive 
the stellar mass function of central galaxies using the condi- 
tional luminosity function (CLF), which expresses how many 
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galaxies of luminosity L reside, on average, in a halo of mass 
M. Using both the galaxy luminosity function and the lu- 
minosity dependence of the correlation length of the galaxy- 
galaxy correlation function obtained from the 2-dcgrcc Field 
Galaxy Redshift Survey (2dFGRS, Colless et al. 2001), Yang 
ct al. (2003) and van den Bosch et al. (2003b) were able to 
put tight constraints on the CLF (see also Yang et al. 2005 
and van den Bosch et al. 2005). As shown in Yang et 
al. (2003), the CLF also allows one to compute the aver- 
age relation between halo mass and the luminosity of the 
central galaxy (assumed to be the brightest galaxy in the 
halo). We have determined this relation using the fiducial 
CLF model considered in van den Bosch et al. (2005; Model 
6 listed in their Table 1). To obtain a stellar mass func- 
tion for the central galaxies, we convert the 2dFGRS bj- 
band luminosity into a stellar mass using a stellar mass-to- 
light ratio M star /L = 4.0 (L/L*) 0,3 (M/L) Q for L < L* and 
M s t ar /L = 4.0(M / L)q for L > L* , which matches the mean 
relation between stellar mass and blue-band luminosity ob- 
tained by Kauffmann et al. (2003) . The resulting stellar mass 
function of central galaxies is shown in Fig. 1(a) as the short- 
dashed curve. For comparison, we also plot the stellar mass 
functions of all galaxies obtained by Bell et al. (2003) (dotted 
curve) and Panter, Heavens & Jimenez (2004) (long-dashed 
curve). Given the large uncertainties involved (see Bell et 
al. 2003 for a detailed discussion), these stellar mass func- 
tions are in remarkably good agreement with each other, 
particularly for the low mass galaxies that are the focus of 
this study. The good agreement suggests that most low mass 
galaxies are indeed central galaxies in small haloes, i.e. satel- 
lite galaxies do not dominate the stellar mass function (see 
also Cooray & Milosavljevic 2005). 

In addition to the stellar mass function, we also con- 
strain our models using the HI mass function of galaxies. 
With large, blind 21-cm surveys that have recently been 
completed, this Hi-mass function has now been estimated 
quite accurately over a relatively large range of masses (see 
Zwaan et al. 2005 and references therein). In Fig. 1(b), we 
show the recent results obtained by Zwaan et al. (2005) and 
Rosenberg & Schneider (2002). Both HI- mass functions are 
well fit by a Schechter (1976) function with a power-law 
slope at the low mass end of about —1.3 ± 0.1. Note that 
this is slightly steeper than the power-law slope at the low 
mass end of the stellar mass function, which is about —1.16 
(Panter et al. 2004). Since galaxy formation is a process 
that involves both stars and cold gas, a combination of ob- 
servational constraints on the luminosity function and the 
Hi-mass function provides important constraints on star for- 
mation and feedback. In fact, as we show below, the Hi-mass 
function constraints are more generic, and it is only by in- 
cluding them we are able to argue against the standard su- 
pernova feedback model. 

2.2 The standard model 

In the standard picture of galaxy formation (White & Rees 
1978) it is assumed that gas cooling conserves specific angu- 
lar momentum. As a result, the baryons cool to form a cen- 
trifugally supported disk galaxy (Fall & Efstathiou 1980). In 
what follows we investigate the mass functions of the cold 
gas and stars of disk galaxies that form within this picture. 
We make the simplifying assumption that each dark mat- 



ter halo forms a single disk galaxy. Clearly this is a severe 
oversimplification since it is known that haloes, especially 
more massive ones, can contain more than one galaxy and 
not every galaxy is a disk galaxy. However, we are only in- 
terested in the properties of low mass haloes, which to good 
approximation contain a single, dominant disk galaxy. In 
particular, the studies of van den Bosch et al. (2003b), Yang 
et al. (2005), and Weinmann et al. (2005) show that in haloes 
with M < 10 12 ft _1 Mq, the mass range considered here, the 
fraction of late-type galaxies is larger than 60 percent (see 
also Berlind et al. 2003). Thus, our simplified model will 
overpredict the number density of disk galaxies in low mass 
haloes, but not by more than a factor of two. Furthermore, 
we will conclude below that the main problem with the stan- 
dard model is one of gas cooling, a problem that will exist 
independent of galaxy type. 

To model the detailed structure of individual disk galax- 
ies we use the model of Mo, Mao & White (1998, here- 
after MMW), which matches a wide variety of properties 
of disk galaxies. Specifically, this model assumes that (i) the 
baryons have the same specific angular momentum as the 
dark matter, (ii) they conserve their specific angular mo- 
mentum when they cool, (iii) they form an exponential disk, 
and (iv) the halo responds to the gas cooling by adiabatically 
contracting. Assumptions (i) and (iv) are supported by nu- 
merical simulations (van den Bosch et al. 2002; Jesseit, Naab 
& Burkert 2002) , while assumption (ii) is required to obtain 
disks of the right size. Finally, assumption (iii) is equivalent 
to assuming a particular distribution of specific angular mo- 
mentum in the proto-galaxy. Haloes are modelled as NFW 
spheres (Navarro, Frenk & White 1997) with a concentration 
that depends on halo mass following Bullock et al. (2001a), 
and a halo spin parameter, A, that is drawn from a lognor- 
mal distribution with a median of ~ 0.04 and a dispersion of 
~ 0.5 (e.g., Warren et al. 1992; Cole & Lacey 1996; Bullock 
et al. 2001b). The one free parameter in this model is the 
disk mass fraction rrid, defined as the disk mass divided by 
the total virial mass. Since radiative cooling is very efficient 
in low mass haloes with M < 10 12 /i _1 Mq, we start our 
investigation by naively setting mj equal to the universal 
baryon fraction, i.e., rrid — 0.17. 

In a seminal paper, Kennicutt (1989) showed that 
star formation is abruptly suppressed below a critical sur- 
face density. This critical density is close to that given by 
Toomre's stability criterion 

Eorit (fl) = (1) 

where k(R) is the epicyclic frequency, <T gas is the velocity 
dispersion of the cold gas and Qcrit ~ 1 (Toomre 1964). 
This critical density determines the fraction of the gas that 
can form stars (Quirk 1972). Given the surface density of 
the disk, Edisk, obtained using the MMW model described 
above, the radius 7?sf where the density of the disk equals 
S cr it can be calculated. Following van den Bosch (2000) we 
assume that the disk mass inside this radius with surface 
density E d i s k > E cr it turns into stars, i.e. 
rRsF 

M star = 2tt / [E disk (i?) - E crit (-R)] RdR (2) 
Jo 

Kennicutt (1989) shows that cr gas = 6km s~ x and Q C rit ~ 1.5 
yields values of Rsf that correspond roughly to the radii 
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Figure 1. (a) The stellar mass functions predicted by the standard model, the model with heating by the UV background, and the 
preheating model, as labelled. The curve labelled 'SN' is the prediction of a model in which cold gas is heated by supernova explosions (see 
text for details). In addition, we show the observational results of Bell ct al. (2003; dotted curve) and Pantcr et al. (2004; long-dashed 
curve), as well as the stellar mass function of central galaxies in dark matter haloes derived from the CLF as described in the text 
(short-dashed curve), (b) The Hi-mass functions predicted by the same three models compared to the observed HI mass functions of 
Rosenberg & Schneider (2002; RS, open circles) and Zwaan ct al. (2005; solid dots). 



where star formation is truncated. However, Hunter et al 
(1998) show that in low surface brightness galaxies Q cr it ~ 
0.75. Therefore, to be conservative, we adopt er gas = 6km s _1 
and Qcrit ~ 0.75. The assumption that all the gas with 
Edisk > E cr it has formed stars is consistent with both obser- 
vations (Kennicutt 1989; Martin & Kennicutt 2001; Wong 
& Blitz 2002; Zasov & Smirnova 2005) and with predictions 
based on the typical star formation rate in disk galaxies 
(Kennicutt 1998) and their formation times (see van den 
Bosch 2001). We compute the gas mass of each model galaxy 
using Mgas = (Mdisk — M star ) where M d i s k = m d M. In other 
words, we assume that the gas surface density E gas = E cr it 
inside Rsf and E gas = Edisk outside 7?sf- Finally, again 
to be conservative, we assume that the molecular hydrogen 
fraction is 1/2 (e.g. Keres et al 2003; Boselli et al. 2002) so 
that the final HI mass of each galaxy, Mm = 0.71M gas /2 
where the factor of 0.71 takes into account the contribution 
of helium and other heavier elements. 

Using the halo mass function given by the ACDM con- 
cordance cosmology, and assuming that each halo hosts a 
single disk galaxy whose properties follow from M, m d , and 
A as described above, we obtain the HI and stellar mass 
functions shown as the solid curves labeled 'standard' in 
Figure 1. Since we are only interested in low-mass haloes, 
and since our model does not include any processes that 
may affect gas assembly in massive haloes, we artificially 
truncate the halo mass function at M = 5 x W 12 ^ 1 M Q , 
which explains the abrupt turnover of the model's stellar 
mass function at high masses. Not surprisingly, our naive 
model severely overpredicts the stellar mass function, yield- 
ing an abundance of systems with M star — W 8 h~ 2 M that 
is two orders of magnitude too large. In addition, the model 



predicts an HI mass function at Mm < W s h~ 2 Mq that is 
more than 10 times higher than the data. Even if we reduce 
the number density of dark matter haloes by a factor of 2, to 
account for the possibility that some isolated haloes may not 
host disk galaxies, the discrepancy remains a factor of five. 
Note that one might try to lower the stellar masses in our 
model by increasing E cr i t , but this leads to an increase of the 
HI masses, which are already too large. Similarly, a decrease 
of E cr it may improve the fit of the HI mass function, but at 
the expense of worsening the fit to the stellar mass function. 
The failure of the standard model to simultaneously fit the 
HI mass function and the stellar mass function is robust to 
the details of star formation. Fitting both mass functions 
simultaneously requires either a modification of the cosmo- 
logical parameters or additional physics to lower mj. In what 
follows, we consider both these possibilities separately. 

2.3 Cosmological Parameters 

One of the main reasons that the predicted HI mass func- 
tion is very steep at the low-mass end is that the halo mass 
function predicted by the standard ACDM model is also 
very steep at the low-mass end. Hence, one way to alleviate 
the discrepancy between the model and observations is to 
change the cosmological parameters such that the low-mass 
slope of the halo mass function becomes shallower. Unfortu- 
nately, the steep halo mass function is a very generic prop- 
erty of all CDM models. In particular, the slope at the low- 
mass end is almost independent of cosmological parameters, 
including the cosmic density parameter and the amplitude 
of the primordial perturbations. The only way to change 
the slope of the mass function at the low-mass end is to as- 
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sume that the effective power index of the primordial density 
perturbation spectrum is significantly lower than the scale 
invariant value. However, such models are not favored by 
the power spectrum derived directly from the temperature 
fluctuations of the cosmic microwave background and the 
Lyman-Q forest (e.g., Croft et al. 1999; Seljak et al. 2003), 
and are difficult, if not impossible, to reconcile with infla- 
tion. Another possibility is that the universe is dominated 
by warm dark matter (WDM) instead of CDM, so that the 
power spectrum on small scales is suppressed by free stream 
damping of the WDM particles. Here again observations of 
the Lyman-a forest provide a stringent constraint on the 
particle mass (Narayanan et al. 2000). With particle masses 
in the allowed range, the WDM model yields a halo mass 
function that is virtually indistinguishable from that of the 
CDM models in the halo mass range considered here. In 
other words, within the parameter space allowed by the 
data, modifications of the cosmological parameters do not 
have any significant impact on the results of the standard 
model presented above. 

2.4 Supernova Feedback 

Thus far we have only considered models in which the disk 
mass fraction is equal to the universal baryon fraction, i.e. 
where = 0.17. We now consider physical mechanisms 
that can lower mj, and investigate whether this allows a si- 
multaneous match of the HI and stellar mass functions. We 
first consider what has become the standard mechanism, 
namely feedback by supernovae. In this model each halo ac- 
quires a baryonic mass fraction that is equal to the universal 
value (0.17) but mj is reduced since supernovae inject large 
amounts of energy into the cold gas, causing it to be ejected 
from the galaxy. 

In semi-analytical models of galaxy formation, two 
schemes have been proposed to model this supernova feed- 
back. In the 'retention' model considered by Kauffmann et 
al. (1999), Cole et al. (2000), Springel et al. (2001) and 
Kang et al. (2005), among others, part of the cold gas in 
a galaxy disk is assumed to be heated by supernovae to the 
halo virial temperature and is added to the hot halo gas for 
cooling in the future. Since radiative cooling is effective in 
galaxy haloes, the feedback efficiency must be high enough 
to keep a large amount of the gas in the hot phase (Benson 
et al. 2003) . However, if there is a critical surface density be- 
low which star formation ceases, no galaxy disk is expected 
to have a surface density below this critical density, because 
otherwise the feedback from star formation would be insuf- 
ficient to keep most of the gas hot. We plot an example of 
such a model in Figure 1. In this model we make the stan- 
dard assumption that the rate at which cold gas is heated 
by supernova explosions is proportional to the star forma- 
tion rate, M rch eat = where /3 = (Vhot/V c ) 2 , with V c 
the circular velocity of the host halo, and Vhot an ajustable 
parameter (e.g. Benson et al. 2003) . If the heated gas is not 
able to cool, the mass of the gas that can form stars will 
be 1/(1-1-/3) times iridM — M gas , where M gas is the mass of 
the gas that remains in the disk. The curve labelled 'SN' in 
Figure 1 is the result corresponding to Vhot = 280kms _1 . 
Although the stellar mass function is reduced significantly 
in this model, the predicted slope at the low-mass end is 
much too steep. Furthermore, the Hi-mass function is un- 



changed from the standard no-feedback model and so it is 
still unable to match the observed Hi-mass function. In ad- 
dition, this model predicts fairly extensive haloes of hot, 
X-ray emitting gas, which is inconsistent with observations 
(Benson et al. 2000). 

An alternative feedback model, considered by Kauff- 
mann et al. (1999) and Somerville & Primack (1999), is the 
'ejection' model in which the reheated gas is assumed to be 
ejected from the current host halo. If the initial velocity of 
the ejected gas is not much larger than the escape velocity 
of the host halo, the gas will be recaptured at a later epoch 
as the halo grows more massive by accreting new material 
from its surroundings. Then , as in the retention model, the 
baryon fraction in the more massive halo will be similar to 
the universal value, except that a (typically small) delay 
time is added. Consequently, this model also cannot pro- 
duce disks with gas surface densities below the critical den- 
sity, which again results in a severe overproduction of low HI 
mass systems. Only if the supernova explosion energy heats 
the cold gas to an energy much greater than the binding 
energy of the halo can the wind escape the halo forever and 
potentially reduce the number of low HI mass systems. How- 
ever, there are several problems with this scenario. First, as 
shown by Martin (1999) and Heckman et al. (2000), the 
observed mass outflow rate in starburst galaxies, which are 
extreme systems, is about twice the star formation rate. This 
implies that one can never achieve a disk mass fraction that 
is lower than about 1/3 the universal baryon fraction, which 
is not nearly enough to match the HI observations. Second, 
the numerical simulations of Mac- Low & Ferrara (1999) and 
Strickland & Stevens (2000) have shown that the mass loss 
rates in quiescent disk galaxies are much lower than those 
observed in starburst galaxies. Third, as shown in Benson et 
al. (2003) , the feedback efficiencies that are required to per- 
manently eject the gas are completely unphysical. Finally, as 
shown in van den Bosch (2002) , even if one ignores all these 
problems and simply ejects the gas forever, the presence of 
a star formation threshold density still assures that the final 
surface density of the gas is similar to the critical surface 
density. As is evident from Figure 5 in that paper, super- 
nova feedback that is modelled with permanent ejection has 
a drastic impact on the stellar masses but leaves the gas 
mass basically unchanged. Again, this owes to the fact that 
gas ejection requires supernovae, which in turn requires star 
formation, which requires a gas surface density that exceeds 
the critical density. Note that although supernova feedback 
may temporarily deplete the gas surface density below the 
critical value, the ongoing cooling of new and previously ex- 
pelled material will continue to increase the cold gas surface 
density until it exceeds E cr i t , initiating a new episode of star 
formation and its associated feedback. As shown in van den 
Bosch (2002), this results in a population of disk galaxies 
whose cold gas surface densities are, in a statistical sense, 
similar to £ C rit- 

In summary, although supernova feedback may be 
tuned to yield a good match to the low-mass end of the 
stellar mass function, it has three fundamental problems: 
(i) the efficiency needed seems unphysically high compared 
to what it achieved in detailed numerical simulations, (ii) 
unless the hot gas is expelled from the halo indefinitely it 
predicts haloes of hot, X-ray emitting gas which are incon- 
sistent with observations, and (iii) as demonstrated here, 
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if one takes account of a star formation threshold density, 
which is strongly supported by both theory and observa- 
tions, it overpredicts the abundance of systems with low HI 
masses by almost an order of magnitude. In short, the prob- 
lem of matching the observed HI mass functions is one of 
preventing gas from entering the galaxy in the first place. 
Standard feedback schemes fail because even if all the gas is 
temporarily removed, e.g. by a massive supernova outflow, 
it will just reaccrete more gas. This is also why our argu- 
ments, although framed around disk galaxies, should hold 
for all galaxies. It is, therefore, important to seek other so- 
lutions that are physically more plausible. 

2.5 Reionization and Preheating 

In the models considered above we assume that the IGM 
accreted by the dark matter haloes is cold, allowing all of 
the gas originally associated with the halo to be accreted 
eventually. However, if the gas in the IGM is preheated to 
a specific entropy that is comparable to or larger than that 
generated by the accretion shocks associated with the for- 
mation of the haloes, not all of this gas will be accreted 
into the halo (e.g. Mo & Mao 2002; Oh & Benson 2003; 
van den Bosch et al. 2003a). In that case, some disks may 
start with a gas surface density already below the critical 
density, making their HI gas masses smaller. Note that this 
circumvents the problem with the supernova feedback sce- 
nario that requires star formation and its inherent high gas 
surface densities. 

Let us first consider photoionization heating by the UV 
background. After reionization, the UV background can heat 
the IGM to a temperature of roughly 20, 000K. As first 
pointed out by Blumenthal et al. (1984), such heating of 
the IGM can affect gas accretion into dark matter haloes of 
low masses (see also Rees 1986; Efstathiou 1992; Quinn et 
all 1996; Thoul & Weinberg 1996). Recent simulations (e.g., 
Gnedin 2000; Hoeft et al. 2005) follow the detailed evolution 
of the UV background and show that at the present time the 
fraction of gas that can be accreted into a dark matter halo 
of mass M can be written approximately as 

mgas = (1 + M C /MY ' (3) 

where fs is the universal baryon fraction. Following Hoeft et 
al. (2005), we consider a model with M c = 1.7 x 10 9 /i _1 M Q , 
and a = 3. Using the same disk formation model as de- 
scribed in Section 2.2, and assuming that disks with surface 
densities below S cr it do not form any stars, we can predict 
the cold gas and stellar masses. The resulting stellar and 
Hi-mass functions are shown as the solid lines labelled 'UV 
in Figure 1. Although preheating by the UV background 
clearly reduces both the stellar and Hi-mass functions at 
the low-mass end, the strength of the effect fails to reconcile 
the models with the data (see also Benson et al. 2002). 

However, since the predicted cold gas mass is signif- 
icantly reduced in haloes with masses below the charac- 
teristic mass scale M c , this motivated us to consider a 
model in which the IGM around low-mass haloes is pre- 
heated to a temperature that is significantly higher than 
20, 000K. A higher temperature corresponds to a larger M c . 
As shown in Lu & Mo (2005, in preparation), the mass frac- 
tion of baryons that are accreted by dark matter haloes 




log[M star ](h-2M ) 

Figure 2. The cold gas mass fraction, defined as the ratio be- 
tween the mass of cold gas and the total mass of stars and cold 
gas, as a function of stellar mass predicted by the preheating 
model (crosses). The thick solid line shows the observed mean 
relation given by McGaugh & dc Blok (1997). 

in a strongly preheated medium is well described by equa- 
tion (3) with q ~ 1. Therefore, as a test of this idea we 
adopt a — 1 and choose M c = 5x 10 n /i _1 M , which corre- 
sponds to an initial specific entropy for the preheated IGM 
of s = T/rie^ 3 ~ lOkevcm 2 (where n e is the number density 
of electrons) . To determine the relationship between M c and 
s we assume that T is the virial temperature of a halo with 
mass M c at the present time and that n e is the mean density 
of electrons within the halo. In §3, we propose a model that 
explains how the IGM is preheated to such a level. Here 
we examine how this preheating affects both the Hi-mass 
function and the stellar mass function 

The solid curves labelled 'preheat' in Figure 1 show the 
stellar and Hi-mass functions predicted by this model, using 
the same disk formation model described above. Contrary to 
the standard model and the reionization model, this preheat- 
ing model agrees with the data fairly well for the low mass 
haloes that concern us here. We underpredict the HI mass 
function at higher masses but remember to be conservative 
we tried to make the HI mass function as small as possible. 
For example, we took a small value Q cr it = 0.75, which is ap- 
propriate for dwarf galaxies. If we took Qcrit ~ 1-5, which is 
appropriate for larger galaxies, it would increase the masses 
and make them better match the observations. 

Unlike the supernova feedback model, preheating can si- 
multaneously match the HI and stellar mass functions. How- 
ever, this does not guarantee that the model also predicts 
the correct ratio of cold gas mass to stellar mass in individ- 
ual galaxies. To test this, we compute for each model galaxy 
the cold gas mass fraction, / gas = M gas /(M gas + M sta r). Fig- 
ure 2 plots / gas as a function of the stellar mass. The scatter 
in the model predictions results from the scatter in the halo 
spin parameters, and is comparable to the observed scatter 
(McGaugh & de Blok 1997). The preheating model predicts 
that the gas mass fraction decreases with stellar mass, in 
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qualitative agreement with the observations (McGaugh & 
de Blok 1997; Garnett 2002). As already demonstrated in 
van den Bosch (2000), this success is a direct consequence 
of implementing a critical surface density for star formation. 
Note that the model predicts gas fractions that are slightly 
lower than those observed. However, given the uncertainties 
involved, both in the data and in the model, and given the 
relatively large amount of scatter, we do not consider this a 
serious shortcoming. 

As a final test of the preheating model we consider 
gas metallicities. The higher gas mass fractions in smaller 
haloes implies that the metallicity of the cold gas must be 
lower in lower mass systems, which is consistent with ob- 
servations (Garnett 2002). However, observations also show 
that the effective metal yield decreases with galaxy lumi- 
nosity (Garnett 2002; Tremonti et al. 2004), suggesting that 
some metals generated by stars must have been ejected 
from the galaxies and that the ejected fraction is larger for 
fainter galaxies. This may seem problematic for the preheat- 
ing model considered here, since we require no outflows to 
match the HI and the stellar mass functions. However, this 
does not exclude the possibility that significant amounts of 
metals have been ejected from low-mass galaxies. In fact, 
as the numerical simulations of Mac- Low & Ferrara (1999) 
demonstrate, supernova feedback in quiescent disk galaxies 
is far more efficient at ejecting metals than mass. The rea- 
son for this is that the metals are largely produced by the 
supernovae themselves, so they are part of the hot bubbles 
of tenuous gas that make up the galactic winds. When these 
bubbles rupture owing to Raleigh- Taylor instabilities this 
strongly metal enriched material, which has relatively little 
mass, is blown away from the disk by its own pressure. Thus, 
although clearly more work is needed to test this in detail, 
we believe that the observed effective metal yields are not 
at odds with the preheating model. 



3 PREHEATING BY GRAVITATIONAL 
PANCAKING 

In the previous section we have shown that a preheating 
model, in which the gas surrounding present-day low mass 
haloes is preheated to a specific entropy of s ~ lOkev cm 2 , 
can simultaneously match the low mass ends of both the HI- 
mass function and the stellar mass function. Here we propose 
a physical process that can cause such preheating. 

We base our proposed model on the following consider- 
ations. In the basic picture of structure formation in CDM 
cosmogonies, as the universe expands, larger and larger ob- 
jects collapse owing to gravitational instability. The col- 
lapse is generically aspherical (Zel'dovich 1970), first form- 
ing sheet-like pancakes (first axis collapse), followed by fil- 
amentary structures (second axis collapse), and eventually 
virialized dark matter haloes (third axis collapse). Thus, ac- 
cording to the ellipsoidal collapse model, the formation of 
a virialized halo requires the collapse of all three axes (e.g. 
Bond & Myers 1996; Sheth, Mo & White 2001). However, 
owing to the large scale tidal field, the density threshold 
for the formation of a low-mass halo, i.e. one with a mass 
below the characteristic mass, M», defined as the mass at 
which the rms fluctuation is equal to unity, can be much 
higher than that in the spherical collapse model. This de- 



lays the formation of low mass haloes relative to the pre- 
diction of the spherical collapse model. Conversely, the tidal 
field accelerates the collapse of the first (shortest) axis and 
hence the formation of a pancake can require a density 
threshold much lower than that for spherical collapse. Con- 
sequently, many of today's low mass haloes, i.e. those with 
M < M,(z = 0) ~ 10 13 /i _1 M , formed in pancakes of 
larger mass, which formed before the haloes themselves. 

In the process of pancake formation, the gas associated 
with the pancake is shock heated. If the temperature of the 
shocked gas is sufficiently high, and if the gas is not able to 
cool in a Hubble time, the haloes embedded in the pancakes 
will have to accrete their gas from a preheated medium, 
which as we showed in the previous section may have impor- 
tant implications for the formation of galaxies within those 
haloes. To see if this process of "gravitational pancaking" 
(or previrialization, Peebles 1990) can generate a preheated 
medium with the required specific entropy, we need to exam- 
ine the properties of the pancakes within which present-day 
low mass haloes formed, and to understand how the gas as- 
sociated with these pancakes was shock heated. 

To study the first problem it is important to realize 
that the bias parameters of haloes with M<0.1M„ have 
similar values (Mo & White 1996; 2002; Jing 1999; Sheth & 
Tormen 1999; Sheth, Mo & Tormen 2001; Seljak & Warren 
2004; Tinker et al. 2004). This means that all such haloes 
are, in a statistical sense, embedded within similar large- 
scale environments. At z = 0, M* ~ lO 13 ^^ 1 Mg, and so all 
haloes with M<Vd 12 h' x M are embedded within similar 
environments. 

According to the ellipsoidal collapse model, the collapse 
of a region on some mass scale M in the cosmic density field 
is specified by 5, the average overdensity of the region in 
consideration, and by e and p, which express the ellipticity 
and prolateness of the tidal shear field in the neighborhood 
of that region. According to Sheth et al. (2001), the density 
threshold for the formation of a virialized halo is given by 
solving 



S cc (e,p) 



1 + /3 



A 2 

5(e ±p ) — 



(4) 



where f3 — 0.47, 7 = 0.615, and S sc is the critical overdensity 
for spherical collapse. For a Gaussian density field, the joint 
distribution of e and p for given a S is 



1 m 1125 < 2 
g(e,p\o) = -j=e (e 



-55 2 (3e 2 +p 2 )/2<T 2 



(5) 



where a is the rms fluctuation of the density field on the 
mass scale in consideration (Doroshkevich 1970). For all e, 
this distribution peaks at p = 0, and the maximum occurs 
at 

e ma x(p = 0\S) = . (6) 

Thus, the most probable value of e is related to the mass 
scale through a. Using this relation, one can obtain a relation 
between the density threshold for collapse and the halo mass: 



6 ec (M,z) = S sc (z) { 1 + 0.47 



(Sheth et al. 2001). 
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Figure 3. The mass [panel (a)], overdensity [panel (b)], gas temperature [panel (c)], and gas specific entropy [panel (d)] of pancakes 
around low-mass haloes at the time of first axis collapse. Thick solid curves assume that perturbations around low-mass haloes have e 
values equal to the most probable values corresponding to the mass scale in consideration, while thick dashed curves show the results in 
which e is assumed to be a constant, i.e. independent of pancake mass. The thin lines in panel (c) show the loci of t coo l = *H for the 
three indicated values of the overdensity. 



The ellipsoidal collapse model can also be used to de- 
termine the density threshold for the formation of pancakes. 
Based on similar considerations, one can obtain a corre- 
sponding relation between the collapse density threshold and 
the mass of the pancake: 



6 ec {M,z) =S sz (z) { 1-0.56 



(8) 



(Shen, Abel, Mo & Sheth, in preparation). As one can see, 
for low peaks (i.e. 8 sc /a <JC 1), the two thresholds can be 
very different, while for high peaks they are similar. Thus, 
the effect of pancaking on subsequent halo formation is more 



important for lower peaks, i.e. for lower mass haloes with 
later formation times. 

Given the above properties of the collapse thresholds, 
we are able to address the following question: for low-mass 
haloes identified at the present time, what is the nature of 
the pancakes within which their progenitor haloes were em- 
bedded at an earlier time? To quantify the formation of pan- 
cakes around a given halo rigorously, one needs to calculate 
the conditional probability distribution for the overdensity 
of density perturbations on various scales around the halo 
and the corresponding tidal shear fields. It is beyond the 
scope of this paper to present such a detailed analysis here. 
Instead, we use the cross-correlation between haloes and the 
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linear density field to determine the average linear overden- 
sity around dark matter haloes on different scales. We then 
use this overdensity to characterize the mean environment 
of haloes of a given mass at the present time at earlier times. 

According to the halo bias model (Mo & White 1996), 
the average linear overdensity within a radius r of a halo of 
mass M can be written as 



6 l (r) = b(M)$ m (r), 



(9) 



where b(M) is the bias parameter for haloes of mass M, 
and the average two-point correlation function of the linear 
density field 



£,m(r')r' 2 dr' , 



(10) 



where £ m (r) is the two-point correlation function. The ra- 
dius r corresponds to a mass scale M r — (47rr 3 /3)p , where 
p is the mean density of the universe at z — 0. As we men- 
tioned above, the bias parameter for present day haloes with 
M < 10 12 /i _1 M is independent of halo mass, with b ~ 0.65 
(Sheth et al. 2001; Jing 1999; Seljak & Warran 2004; Tinker 
et al. 2004). We adopt this value to calculate Si. It is then 
straightforward to calculate the linear overdensity Si as a 
function of r and the corresponding mass scale M r . If this 
overdensity reaches the overdensity threshold for pancake 
formation at a given redshift, a pancake of mass M r will 
form. Since the overdensity threshold depends on the values 
of p and e, in principle one has to calculate the joint dis- 
tribution function of c and p for the appropriate mass and 
overdensity, under the condition that the region contains a 
halo of mass M at the present time. As an approximation, we 
assume that both e and p take their most probable values on 
the mass scale in question, so that p = and e = e max . We 
expect this approximation to be valid for M r 3> M, where 
the correlation between the properties of the halo and its 
environment becomes weak (e.g. Bardeen et al. 1986). 

For a given 8i, M r , p and e, we use the ellipsoidal col- 
lapse model described in Bond & Myers (1996) to follow the 
collapse of a uniform ellipsoid embedded in an expanding 
background along all three axes, taking into account both 
the ellipsoid's self-gravity and the external tidal field. Fol- 
lowing Bond & Myers, we assume that each axis freezes out 
at a constant fraction of its initial radius, so that the mean 
overdensity of the collapsed object is just the same as that 
in the spherical collapse model (see Bond & Myers 1996 for 
details). The solid curve in Fig. 3(a) shows the mass of the 
pancake that forms around present day low-mass haloes, as 
a function of z. Remember that at earlier times these haloes 
have smaller masses. Also, owing to the roughly constant 
bias parameter, these results are almost independent of halo 
mass for M < 1O 12 /i _1 M0. As one can see, the pancake mass 
decreases with redshift, because the overdensity threshold 
for collapse is higher at higher z. At z ~ 2, the pancake 
mass is about 10 12,5 /i _1 M . The solid curve in Fig. 3(b) 
shows the overdensity of the pancake at the time of forma- 
tion. This overdensity increases with z, and is about 10 for 
z = 1-2. The ellipsoidal collapse model also determines the 
velocity along the first axis at the time of pancake forma- 
tion. The gas associated with the pancake will be shocked. If 
we assume an adiabatic equation of state and that the shock 
is strong, all the kinetic energy is transformed into internal 
energy. We can calculate the post shock gas temperature as 



(11) 



16fc s ' 

where ks is Boltzmann's constant, n is the mean molecular 
weight, and V\ is the velocity along the first axis at the time 
of shell crossing. We plot this temperature, as a function 
of 2, as the thick solid curve in Fig. 3(c). The temperature 
decreases with increasing redshift because the mass of the 

1/3 

pancake, M p , is smaller at higher z and V\ ~ HqR p oc M p 
where R p is the Lagrangian radius of the pancake. At z ~ 2, 
the temperature is about 10 5,5 K. Assuming that the gas 
overdensity is the same as that of the dark matter, we can 
estimate the specific entropy generated in the shock as 



T 

2/3 



= 17 X 



ft B , h 2 



-2/3 



0.024 



(io 55 k) 



keV cm 



(12) 



where we have taken /i — 0.6, valid for a completely ion- 
ized medium dominated by hydrogen and helium. As shown 
by the solid curve in Fig. 3(d), s increases rapidly with de- 
creasing redshift, mainly owing to the decreasing gas den- 
sity. At z ~ 2 the resulting entropy is s ~ 15KeV cm 2 . For 
z < 2.5, pancake formation results in a preheated IGM with 
s > lOKeV cm 2 , which corresponds to the value adopted in 
the preheating model discussed in the previous section. 

The results presented here are based on the assumption 
that both e and p (which specify the local tidal field) take 
their most probable values. In reality, the tidal field around 
a point in a Gaussian density field must be coherent over 
a finite scale. Since the low-mass haloes at z = are low 
peaks, typically with large values of e, it is possible that 
the value of e for the region around such a halo is larger 
than the most probable value. Without going into detailed 
calculations that include correlations of the tidal field on 
different scales, we consider an extreme case in which we 
set e = 0.45 for all M r . This value for e is approximately 
equal to that for a peak with S/a = 1. The thick dashed 
curves in the four panels of Fig. 3 show the results for this 
extreme model. Because the assumed value of e is larger 
than the most probable value, a given mass pancake forms 
earlier and correspondingly the overdensity of the pancake 
is lower. However, for z < 3, the change in T is less than 50% 
and the change of s is less than a factor of two. 

The preheating by gravitational pancaking is expected 
to have important consequences for the subsequent accretion 
of gas into dark matter haloes only if the heated gas cannot 
cool efficiently. The cooling time of the heated gas can be 
written as 



6.3 Ali, 



Mb, oh 2 
0.024 



—) 

10 55 K/ 



(13) 



where A_23 is 

the cooling function in units of 10 -23 erg sec -1 cm 3 . This 
time should be compared with the Hubble time, 2 



2 Our notation is such that H(z) = 1 for an Einstcin-de Sitter 
Universe, which is also valid to good approximation for the ACDM 



concordance cosmology at z>l. 
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where H(z) = fij /2 (1 + zf /2 H /H(z) and H(z) is Hubble's 
constant at redshift z. In Fig. 3(c), the three thin curves 
show the loci of t cooi — £h in the T-z plane for 5 = 5, 10, 
and 20, respectively. Here we have used the cooling function 
of Sutherland & Dopita (1993) with a metallicity of O.OIZq. 
Effective radiative cooling only occurs below these curves. 
Comparing these curves with those showing the temperature 
of the gas in pancakes, we see that heating by gravitational 
pancaking at z < 2 can have a significant impact on the sub- 
sequent accretion of gas into haloes that will be low mass 
today. At higher redshifts, however, the cooling proceeds suf- 
ficiently fast so that the IGM within the pancake can cool 
back to its original temperature by the time the low mass 
haloes within the pancake form. 

Once again remember that the specific entropy gener- 
ated by gravitational pancaking at z ~ 2 is very similar to 
what one needs to suppress the cold gas fraction in low-mass 
haloes (see §2). Thus we conclude that the preheating model 
discussed in the previous section, which is extremely success- 
ful at explaining both the stellar and HI mass functions, has 
a natural origin. One does not need to invoke any star for- 
mation or AGN activity; rather, the IGM is preheated to the 
required specific entropy by the same process that explains 
the formation of the dark matter haloes themselves. 



4 DISCUSSION 

4.1 Comparison with numerical simulations 

In the previous sections we have argued that preheating by 
gravitational pancaking causes the shapes of the HI and stel- 
lar mass functions at the low mass ends to agree with obser- 
vations. This begs the question as to why this effect has not 
already been seen in existing hydrodynamical, cosmological 
simulations. The short answer is that no previous simulation 
had the necessary resolution. 

To study the effects discussed above places two differ- 
ent resolution constraints on simulations. First, they must be 
able to resolve the shocks that occur in the forming pancakes 
and second, they must resolve the small galaxies that form 
within them. The typical pancake thickness is ~ 200h~ 1 kpc 
in comoving units, with an overdensity of about 10. In a 
Smoothed Particle Hydrodynamics (SPH) code, it requires 
4h s to resolve a shock (Hernquist & Katz 1989), where h B 
is the SPH smoothing radius usually chosen so that there 
are 32 particles within a sphere of radius 2h s . Since the 
pancake has a shock on both sides, the absolute minimum 
resolution has to be at least 8h s across the pancake thick- 
ness, i.e. h s < 25h~ kpc. Even a shock capturing Eulerian 
grid or AMR code requires at least 4 grid cells across the 
pancake width to resolve the post-shock structure. The typ- 
ical galaxy that needs to be affected has a dark halo with 
a circular velocity of ~ 50 kms -1 , which corresponds to a 
halo mass of ~ 10 10 Mq and a virial radius of ~ 50/i _1 kpc. 
To marginally resolve such haloes requires at least 100 dark 
matter particles and hence a particle mass for the dark mat- 
ter of less than 10 s Mq . In addition, one must have at least 
2 spatial resolution elements within the virial radius and 
hence the spatial resolution must be at least 25/i _1 kpc. In 



an Eulerian code both these spatial resolution requirements 
are identical but in a Lagrangian code like SPH, where the 
spatial resolution is variable, the halo requirement is actu- 
ally (200/10) 1 / 3 = 2.7 times easier to satisfy owing to the 
higher overdensity. 

The above resolutions have not been achieved in any 
published numerical simulation. The highest resolution sim- 
ulation to date at z = (Keres et al. 2005) has 128 3 gas 
and dark matter particles in a periodic, cubic volume of 
22.22ft -1 comoving Mpc on a side. At an overdensity of 10 
this simulation has h s — 81/i -1 kpc, more than 3 times too 
large to properly resolve the pancake thickness. In addition, 
the dark matter particle mass of 9 x 10 8 Mq is almost an 
order of magnitude too large. Furthermore, the volume is 
really too small to evolve the simulation down to z — 
and to sample enough different environments. Springel & 
Hernquist (2003) have a SPH simulation in a large enough 
volume, 100ft -1 Mpc on a side with 324 3 gas and dark mat- 
ter particles. However, both the spatial and mass resolution 
are worse than in the case of Keres et al. (2005), 143/i _1 
kpc and 2 x 10 9 Mq, respectively. Their latest unpublished 
simulation has 486 3 particles but still has worse spatial reso- 
lution (96/i _1 kpc) than Keres et al. (2005). In addition, the 
particle mass is 6 x 10 Mq , which is still six times too large. 
The Eulerian simulation of Kang et al. (2005b) has 1024 3 
cells in the same size volume as Springel & Hernquist (2003) 
for a grid cell size of 97/i _1 kpc, again worse than Keres et 
al. (2005). Finally, Nagamine et al. (2001) have 768 3 cells 
within a volume of 25h~ comoving Mpc on a side with a 
spatial resolution of 75/i _1 kpc. 

To perform a cosmological SPH simulation in a uniform 
periodic volume large enough to evolve robustly to z = 0, i.e. 
lOO/i -1 Mpc on a side, and marginally resolve preheating, 
as outlined above, would require 1860 dark matter and gas 
particles. Such a large simulation is well beyond the reach 
of the current generation of computers and codes. However, 
one can reach such high resolutions in a large volume by 
using a "zoom-in" strategy, in which the resolution is only 
high in the region of interest (Katz & White 1993). Here, 
one starts with a large volume simulated at moderate res- 
olution, identify the object whose formation one wishes to 
study, trace the particles that end up in or near this ob- 
ject back to the initial conditions, replace the particles in 
this Lagrangian region with a finer grid of less massive par- 
ticles, add the small scale waves that can be resolved by 
this finer grid to the initial fluctuation spectrum, and re-run 
the simulation. The particle density away from the region of 
interest is reduced by sparse sampling the original particle 
grid in a series of nested zones, always keeping the particle 
density high enough to maintain an accurate representation 
of tidal forces. This approach is similar to the one used by 
other groups (e.g. Navarro & Steinmetz 1997, 2000; Stein- 
metz & Navarro 1999, Navarro et al. 1995; Sommer-Larsen 
et al. 1999, Robertson et al. 2004, Governato et al. 2004) 
in their simulations of individual galaxies but one needs to 
focus on present day low mass galaxies that form within 
pancakes at higher redshift. To test our predictions regard- 
ing the preheating by gravitational pancaking we plan to 
carry out investigations along this line in the near future. 
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4.2 Implications for galaxy formation and the 
IGM 

According to the results presented above, the assembly of 
gas into galaxy-sized haloes is expected to proceed in two 
different modes with a transition at z ~ 2. At z > 2, the 
accreted intergalactic gas is cold. Since radiative cooling is 
efficient in galaxy haloes, gas assembly into galaxies is ex- 
pected to be rapid and to be dominated by clumps of cold 
gas. Combined with the fact that the formation of galaxy- 
sized dark matter haloes at z > 2 is dominated by major 
mergers (e.g. Li et al. 2005), this suggests that during this 
period gas can collapse into haloes quickly to form starbursts 
and perhaps also feed active galactic nuclei (e.g. Baugh et 
al. 2005). Galactic feedback associated with such systems 
may drive strong winds into the IGM, contaminating the 
IGM with metals. Note that strong feedback in this phase 
is required to reduce the star formation efficiency in high- 2 
galaxies; otherwise too much gas would already turn into 
stars at high z. This phase of star formation maybe what is 
observed as Lyman-break galaxies and sub-mm sources at 
z > 2, and maybe responsible for the formation of elliptical 
galaxies and the bulges of spiral galaxies. At z < 2, however, 
the situation is quite different. Since the medium in which 
galaxy-sized haloes form is already heated by previrializa- 
tion and since radiative cooling is no longer efficient, the 
accretion is expected to be dominated by hot, diffuse gas. 
Such gentle accretion of gas might be favorable for the for- 
mation of the quiescent disks of spiral galaxies. Because the 
accreted gas is diffuse rather than in cold clumps, it can 
better retain its angular momentum as it settles into a rota- 
tionally supported disk. In addition, since the baryonic mass 
fraction that forms the disk is significantly smaller than the 
universal baryon fraction, the disk is less likely to become vi- 
olently unstable. Both effects may help alleviate the angular 
momentum problem found in some numerical simulations, 
i.e. disks that form in CDM haloes have too little angular 
momentum and are too concentrated (Navarro & Steinmetz 
1997). 

Depending on the halo formation history, the bulge-to- 
disk ratio will vary from system to system. For haloes that 
have assembled large amounts of mass before preheating, 
the galaxies that form within them should contain signifi- 
cant bulges, while in haloes that form after preheating the 
galaxies should be dominated by a disk component. Since 
haloes that form later are less concentrated, this model nat- 
urally explains why later-type galaxies usually have more 
slowly rising rotation curves. An extreme example would be 
low-surface brightness (LSB) galaxies. By definition, LSB 
galaxies are disks in which the star formation efficiency is 
low. These galaxies are also gas rich, have high specific an- 
gular momenta, and show slowly rising rotation curves. The 
last two properties are best explained if LSBs are hosted 
by haloes that have formed only in the recent past, since 
such haloes have low concentrations (e.g. Bullock et al. 2001; 
Wechsler et al. 2002; Zhao et al. 2003a, b), and high spin pa- 
rameters (e.g. Mailer, Dekel & Somerville 2002; D'Onghia & 
Burkert 2004; Hetznecker & Burkert 2005). However, there 
is one problem with such a link between LSBs and newly 
formed dark matter haloes. Since the formation of such 
haloes involves major mergers, these systems are expected 
to produce strong starbursts rather than low-surface bright- 



ness disks. This problem does not exist in our model since 
these haloes are expected to form in a preheated medium 
and gas accretion into such haloes will be smooth. It will be 
interesting to see if our model can predict the right num- 
ber of LSB galaxies and explain the existence of extreme 
systems like Malin I. 

Our model also opens a new avenue to understand the 
evolution of the galaxy population with redshift. As we ar- 
gued above, star formation at z > 2 is expected to be dom- 
inated by starbursts associated with major mergers of gas 
rich systems, while star formation at z < 2 is expected to 
occur mostly in quiescent disks. This has important impli- 
cations for understanding the star formation history of the 
universe and for understanding the evolution of the galaxy 
population in general. For example, our model implies a 
characteristic redshift, z ~ 2, where both the star forma- 
tion history and galaxy population make a transition from a 
starburst-dominated mode to a more quiescent mode. Fur- 
thermore, if AGNs are driven by gas-rich major mergers, a 
transition at z ~ 2 is also expected in the AGN population. 
There are some hints about such transitions in the observed 
star formation history (e.g. Blain et al. 1999), and in the 
observed number density of AGNs (e.g. Shaver et al. 1996). 
Recent observations of damped Lyman alpha systems also 
suggests a change in behavior at z ~ 2 both in the cold gas 
content and in the number density of such systems (see Rao 
2005). 

The preheated medium we envision here is closely re- 
lated to the warm-hot intergalactic medium (WHIGM) un- 
der intensive study in recent years. Hydrodynamical simula- 
tions show that between 30 and 40 percent of all baryons re- 
side in this WHIGM, which is produced by shocks associated 
with gravitational collapse of pancakes and filaments (e.g. 
Cen & Ostriker 1999; Dave et al. 2001; Kang et al. 2005b). 
These results are consistent with observational estimates 
based on the study of UV absorption lines in the spectra of 
low-redshifts QSOs (see Tripp et al. 2004 for a review). In 
our model, the low-temperature component of this WHIGM, 
i.e, with temperatures of a few times 10 5 K, is associated with 
pancakes of relatively low mass (~5x lO 12 ^ 1 M ), within 
which late-type galaxies form. As we discussed in §4.1, cur- 
rent simulations are still unable to make accurate predic- 
tions regarding this particular component of the WHIGM. 
Future simulations of higher resolution, however, may shed 
light on the relation between the properties of galaxies and 
those of the IGM in their immediate surroundings. Such a 
relationship may prove pivotal in observational hunts for the 
missing baryons as the spatial distribution of galaxies and 
their properties can serve as guideposts in the search for the 
WHIGM. 



5 CONCLUSIONS 

Understanding the shallow faint-end slope of the galaxy lu- 
minosity function, or equivalently the stellar mass function, 
is a well known problem in galaxy formation. In the standard 
model it is often assumed that supernova feedback keeps 
large fractions of the baryonic material hot, thus suppress- 
ing the amount of star formation. Although the efficiency 
of this process might be tuned so that one fits the faint- 
end slope of the galaxy luminosity function it has a number 
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of problems. First, the required efficiencies are extremely 
high and are inconsistent with detailed numerical simula- 
tions (e.g., Mac-Low & Ferrara f999). Second, unless the 
hot gas is somehow expelled from the dark matter halo for- 
ever, which requires even higher feedback efficiencies, this 
model predicts hot haloes of X-ray emitting gas around disk 
galaxies, which is inconsistent with observations (e.g., Ben- 
son et al. 2000). In this paper we have demonstrated that 
this model has an additional shortcoming. Using recently ob- 
tained HI mass functions we show that the supernova feed- 
back model predicts HI masses that are too high. The reason 
is that supernova feedback requires star formation which in 
turn requires high surface densities of cold gas. The latter 
owes to the existence of a star formation threshold density, 
which has strong support from both theory (e.g., Quirk 1972, 
Silk 2001) and observations (e.g., Kennicutt 1989; Martin & 
Kennicutt 2001). 

We therefore argue that simultaneously matching the 
shallow, low-mass slopes of the stellar and HI mass functions 
requires an alternative mechanism, which does not directly 
depend on star formation. We demonstrate that a mecha- 
nism that can preheat the IGM to a specific gas entropy of 
~ 10 keV cm 2 , can fit both the observed stellar mass function 
as well as the HI mass function. We also show that gravita- 
tional instability of the cosmic density field can be the source 
of this preheating. Our idea is fairly simple: low mass haloes 
form within larger-scale overdensities that have already un- 
dergone collapse along their first axis. This 'pancake' forma- 
tion causes the associated gas to be shock-heated, and pro- 
viding that the gas cooling rate is sufficiently slow, the low 
mass haloes embedded within these pancakes subsequently 
form in a preheated medium. 

Using the ellipsoidal collapse model, we demonstrate 
that the progenitors of present-day haloes with masses 
M < 10 12 /i -1 M Q were embedded in pancakes of masses 
~ 5 x lO 12 ^ -1 Mq at z ~ 2. The formation of such pan- 
cakes can heat the gas associated with them to a temper- 
ature of ~ 5 x I0 5 K and compress it to an overdensity of 
~ 10. This gas has a cooling time longer than the age of 
the Universe and a specific entropy of about I5keVcm 2 ; the 
amount needed to explain the observed stellar and HI mass 
functions. 

Our results demonstrate that heating associated with 
previrialization may also help solve a number of outstanding 
problems in galaxy formation within a CDM universe. How- 
ever, detailed, high-resolution numerical simulations will be 
required to test our predictions in detail. Such simulations 
will help us understand how the formation of a pancake 
heats the gas initially associated with it and whether the 
amount of heating follows our analytic results. For exam- 
ple, our calculation indicates that shock heating will domi- 
nate over cooling in typical pancakes at z < 2. However, such 
calculations assume that the gas is smoothly distributed. 
Structures at scales smaller than the pancake could lead to 
density inhomogeneities, either pancakes, filaments or ha- 
los, and these could promote extra cooling and move the 
transition to lower redshifts. A similar effect occurs when 
one calculates the transition mass from cooling dominated 
to infall dominated accretion in galaxy formation (White & 
Frenk 1992) and compares it with actual simulations (Keres 
et al 2005). They will also help us understand how such 
heating affects subsequent gas accretion and cooling into 



the dark haloes that form in the pancake and hopefully de- 
rive actual stellar and HI mass functions at the small mass 
end. We have argue that no cosmological, hydrodynamical 
simulation to date has the required spatial and/or mass res- 
olution to study the pancake preheating proposed here. To 
achieve the required numerical resolution, we suggest res- 
imulating, at high resolution, a number of pancakes (and 
filaments) with masses of the order of 5 x lO 12 /!" 1 Mq, iden- 
tified from large cosmological, hydrodynamical simulations. 
Thus far this 'zoom-in' strategy has mainly been used to 
study clusters and galaxies. If our estimates are correct, it 
will be extremely interesting to apply this technique to study 
pancakes and their enclosed, low mass haloes. 
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